First Principles Study of CaFe 2 As 2 "Collapse" Under Pressure 
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We perform first principles calculations on CaFe2 As2 under hydrostatic pressure. Our total energy 
calculations show that though the striped antiferromagnetic (AFM) orthorhombic (OR) phase is 
favored at P=0, a non-magnetic collapsed tetragonal (cT) phase with diminished c-parameter is 
favored for P > 0.36 GPa, in agreement with experiments. Rather than a mechanical instability, 
this is an enthalpically driven transition from the higher volume OR phase to the lower volume cT 
phase. Calculations of electronic density of states reveal pseudogaps in both OR and cT phases, 
though As(p) hybridization with Fe(d) is more pronounced in the OR phase. We provide an estimate 
for the inter-planar magnetic coupling. Phonon entropy considerations provide an interpretation of 
the finite temperature phase boundaries of the cT phase. 

PACS numbers: 74.70.Xa,74.62.Fj,74.20.Pq,74.25.Bt 
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Recent discoveries [l[ of iron-based pnictides provide 
a rich arena to explore the interplay between structural, 
magnetic and superconducting properties, and the conse- 
quent emergence of new physics. These materials provide 
insight into the competing roles of magnetism and pair- 
ing correlations, such as in the high temperature cuprate 
superconductors. Amongst the pnictides, the 122 ternary 
compounds ^4Fe2As2 (A = alkaline earth metal Ca, Ba, 
Sr), belonging to the ThCr2Si2 structure family, draw 
particular interest owing to the rich behavior observed 
upon chemical substitution or applying pressure > 
such as different structural phases and superconductiv- 
ity. Applied pressure has the advantage of introducing 
less disorder compared to chemical substitution. 

CaFe2As2, the smallest- volume member of this family, 
is of great current interest as it serves as a readily ac- 
cessible system that exemplifies the key features of the 
AFe2As2 compounds 0,0- At ambient pressure, at T c i 
= 170 K, CaFe2As2 undergoes a lst-order transition from 
a high temperature tetragonal (T) phase to a low temper- 
ature orthorhombic (OR) phase, that is striped along the 
a-axis and antiferromagnetically (AFM) ordered along 
the c-axis. This may be viewed as a magneto-structural 
transition from a high-T phase with fluctuating mag- 
netic moments [8] to one with long-range AFM order. 
The striped magnetic order drives the orthorhombic sym- 
metry breaking with the antiferromagnetic bonds in the 
a direction being slightly longer than the ferromagnetic 
bonds in the b direction. The T-OR transition tempera- 
ture T c i decreases with applied pressure. 

At low-T, under hydrostatic pressure P ~ 0.35 GPa 0, 
the system undergoes a transition from the AFM-OR 
phase to a non-magnetic tetragonal phase, but with a 
compressed c-axis value; this has been termed the "col- 
lapsed" tetragonal phase (cT). At high-T, and P > 0.35 
GPa, another lst-order transition occurs at T C 2, from 



the tetragonal T to the collapsed cT phase. T C 2 in- 
creases with pressure. Several features are sensitive to 
pressure conditions; in particular, lack of superconduc- 
tivity up to P ~ 0.65 GPa for the case of hydrostatic 
pressure Q, compared to observation of superconduc- 
tivity under conditions creating uniaxial pressure [J, 0] ■ 
Some experiments [nil ] indicate the presence of a low-T 
tetragonal phase sandwiched between the OR and the cT 
phase, suggesting that that superconductivity in a nar- 
row region may be facilitated by the fluctuating moments 
present in the T phase. The transition from the OR to 
the cT phase occurs at lower pressures in the uniaxial 
case 

Prior electronic density functional theory (DFT) 
work [11| has considered the pressure and doping depen- 
dence of BaFc2As2- More recently, DFT studies com- 
pared non-hydrostatic (i.e. anisotropic) and hydrostatic 



pressure dependences of BaFe2As2 and CaFe 2 As2 [12, 
13[, and proposed [l4[ a Hund's rule coupling model of 
the phase transitions in these compounds. Our DFT 
work provides a different understanding of CaFe 2 As2 un- 
der hydrostatic pressure, and goes beyond previous DFT 
work to incorporate a thermodynamic analysis of non- 
zero temperature and pressure. 

The key results of our work are as follows. Total 
energy considerations, as a function of pressure, lead 
to findings that agree well with experimental determi- 
nation of electronic structure and magnetic ordering. 
Though the higher volume AFM-OR phase is favored 
at P = 0, the lower volume non-magnetic cT phase is 
favored for P > 0.36 GPa. Hence we assert the OR- 
cT transition is enthalpy-driven. Our density of states 
(DOS) calculations shows pseudogaps in both OR and 
the cT phases, with the pseudogap being narrower and 
deeper, and As(p) hybridization with Fe(d) more pro- 
nounced, in the OR phase. Considering ferromagnetic 
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FIG. 1: (Color online) Structures used in calculations; chemical species are indicated by color, vertical heights are indicated 
by atomic size, (a) Left panel shows the c-axis view, with afe-plane stripe ordering denoted by +/-; (b) right panel shows the 
6-axis view with the a-axis stripe and the c-axis AFM ordering denoted by +/-. 



(FM) and antiferromagnetic stripe orders, we provide an 
estimate of the inter-planar magnetic coupling. Finally, 
we show that thermodynamic calculations involving lat- 
tice phonons may provide useful insight into the observed 
phase boundaries at non-zero temperature and pressure. 

The structures of interest are tetragonal, with Pear- 
son symbol tllO denoting body centered tetragonal with 
10 atomic sites per unit cell, and orthorhombic, with 
Pearson symbol oF20 denoting face-centered orthorhom- 
bic with 20 atomic sites per unit cell. The oF20 crystal 
structure is a based on a y/2 x \/2 R45 tetragonal super- 
cell of the tllO structure, followed by a weak orthorhom- 
bic distortion. To improve consistency of our calculated 
property differences, we employ this tetragonal super- 
cell for studies of the tetragonal structure, so we include 
20 atoms in all our reported calculations. Calculations 
of the orthorhombic structure utilize spin polarization, 
with initial moments in the striped antiferromagnetic ar- 
rangement. Specifically, spins are ordered antiferromag- 
netically in the a direction, ferromagnetically in the b 
direction (i.e. "striped" in the afo-plane) and antiferro- 
magnetically in the c direction; see Fig. [TJ The symmetry 
group of the atomic positions is IA/mmm (no. 139) for 
tllO, and Fmmm (no. 69) for oF20. 

We utilize VASP 0, O to carry out first principles 
total energy calculations, adopting projector augmented 
wave potentials 17|, LL8J. For a density functional we 
choose the PBE generalized gradient approximation, as 
calculations utilizing LDA fail to stabilize the striped 
antiferromagnetic orthorhombic phase, and the PW91 
GGA predicts the orthorhombic phase to be energeti- 
cally unstable at low temperature and pressure. We relax 
all atomic positions and lattice parameters, and increase 
our k-point densities (to an 8x8x4 Monkhorst mesh) and 
plane-wave energy cutoff (to 340 eV) until energy differ- 
ences have converged to 0.1 meV/atom. Unusual care 
in selection of density functional and convergence is re- 
quired because of the extremely small energy differences 
of order 1 me V/ atom that must be resolved. 



Fig. [5] shows the results of our calculation of total 
energy vs. volume for the nonmagnetic cT and the 
AFM-OR phases, as described above. For each struc- 
ture a we fit E a (V) to a cubic polynomial of the form 
E a = E a ° + (l/2V a )B a (V-V a ) 2 + (l/6V a 2 )C a (V-V a ) 3 , 
where V a and E a ° are the volume and energy at P = 
—dE/dV = 0, B a is the bulk modulus and C a is the 
nonlinear bulk modulus. Fitted values of these quanti- 
ties are listed in Table HI The line of double tangency 
E t = E° - P°V, where P c ° = 0.022 eV/A 3 (0.36 GPa) is 
the critical pressure. 

Fig. [2] implies that though the high volume OR phase is 
energetically favorable at low pressure, the enthalpies (H 
= E + PV) of the phases cross at P®, beyond which the 
high energy cT phase, with diminished c-parameter, has 
lower enthalpy relative to that of the OR phase. This is 
more clearly seen in Fig.[3]where we plot, for the OR and 
cT phases, the enthalpy difference (AH=H c t-Hqr), and 
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FIG. 2: (Color online) Calculated total energies (E) of the 
collapsed tetragonal (cT) and orthorhombic (OR) structures 
plotted against volume (V). The solid lines are cubic polyno- 
mial fits (see text; Table 1); the dotted line of double tangency 
gives the critical pressure beyond which cT is stable despite 
its higher energy. 
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Structure 


b/a c/a Va 




TTlFe 


cT 


1 1.84 16.55 


-6.1296 85.6 -654 





OR 


0.98 2.02 17.57 


-6.1318 63.3 -152 


1.78 



TABLE I: Properties of structures at V=V a . Length and 
volume units are in A, energy in eV/atom, bulk moduli in 
GPa, and magnetic moment in Bohr magnetons. 



the lattice parameters (a, b, c) , versus pressure: While 
both the OR and the cT structures are stable across the 
pressure range studied, AH vanishes at P® = 0.36 GPa, 
causing the c-axis lattice parameter to switch from the 
higher value of 11.45 A of the OR phase to a lower 
value of 10.45 A, the defining feature of the collapsed cT 
phase. Thus, the OR-cT transition at T=0 and finite-P is 
a transition in thermodynamic stability rather than a soft 
mode or irreversible mechanical instability as implied by 
the term "collapse" . In particular, the OR phase contin- 
ues to exist at high pressure, and the cT phase is present 
at low pressure, they are simply metastable rather than 
thcrmodynamically stable. 

Our calculated lattice parameters, shown in Fig. [31 
across the pressure range (0 < P < 1 GPa) are within 
2% of experiments. At P° the discontinuities in lat- 
tice parameters are Aa=+0.04 A, Ab=+0.14 A, and 
Ac=— 1.00 A, in good agreement with experiments 
On extending our calculations to higher pressures, P = 
1.5 GPa, we do not find signatures for any other low-T 
transition. 

In Fig. [4] we show our P = calculated total electronic 
density of states (DOS), as well as the Fe, As, and Ca 
local DOS, for the OR and cT phases. The As(p) hy- 
bridization with Fe(d) is found to be more pronounced 
in the AFM-OR phase than in the non-magnetic cT phase 
as can be seen in the As local DOS within 2 eV below 
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FIG. 3: (Color online) . The top two panels show the behav- 
ior of the calculated lattice parameters, a, b, c of OR and cT 
structures with pressure. The bottom panel shows the change 
in enthalpy, AH vs pressure, P. Note the "jumps" in the lat- 
tice parameters at P® = 0.36 GPa (dotted lines), where AH 
= . 
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FIG. 4: (Color online) Calculated partial and total density 
of states (DOS) in cT (top panel) and OR (bottom panel) 
phases at P = 0. Ef is denoted by pink dashed lines. Note the 
pseudogap features in OR and cT phases. As(p) hybridization 
with Fe(d) can be seen in OR phase. 



Ep , and in the increased Fe d-band width over the same 
range, in the OR phase. Ca(s) orbital participation is 
insignificant in both cases. We have checked the s,p,d 
orbital-decomposed DOS and the band dispersion rela- 
tions (not shown here) for consistency. This overall be- 
havior persists for P>0. Both the OR and cT phases ex- 
hibit pseudogaps. The Fermi level (Ef) lies in the center 
of the pseudogap in the OR phase, while the pseudogap is 
broader and shifted up from Ep in the cT phase. We note 
that the width and location of our DFT-calculated pseu- 
dogaps compare well with LDA+DMFT calculations (i~9j 
at P = 0. 

As in other DFT calculations, our values of the Fe- 
moment at all pressures are large compared to the exper- 
imental results. At P = 0, the calculated mp e ~ 1.78 /is, 
compared to experimental value of 0.8 hb- This decreases 
slightly with pressure. We estimate the inter-planar mag- 
netic coupling based on our calculations of total energy 
and Fe magnetic moment, mp e . The Heisenberg Hamil- 
tonian for nearest-neighbor interactions is 



H = - 



<ij> 



T- ■ S ■ • S ■ 



(1) 



where Jij is the nearest-neighbor Heisenberg coupling 
between neighbors and have opposite signs for fer- 
romagnetic and antiferromagnetic ordering, and Si are 
the spin operators at site i. Under the usual assump- 
tion that the magnetic coupling is site-independent, the 
inter-planar magnetic coupling J± is related to the dif- 
ference in the total energy E of the stripe layers (Fig.QJi) 
stacked ferromagnetically (FM) and antiferromagneti- 
cally (AFM, as shown in Fig. Hb); i.e. AE = E FM - 
Eafm = 2Ji m| e . Our DFT calculations give AE=8.5 
meV/atom. Using the calculated mp e sa 1.78 \ib at P=0, 
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we estimate J± « 2.7 meV. However, magnetization pat- 
terns with tetragonal symmetry (i.e. wnstriped along the 
a-axis, and also checkerboard in the a&-plane) lost their 
magnetization and relaxed to the cT state, illustrating 
the need to supplement eq. ((TJ with coupling of mag- 
netism to lattice relaxation. 

Although the first principles results are derived at T=0 
K, we can estimate the temperature-dependent variation 
of critical pressure P C (T) at low pressure to predict the 
finite-T phase boundary between the orthorhombic and 
collapsed tetragonal phases. The Clausius-Clapeyron 
equation relates the slope of the phase boundary to 
the discontinuities in entropy (S) and volume (V) as 
dP/dT=dAS/dAV. We calculated AV = V cT - V OR « 
1 A 3 /atom. Temperature-dependent entropy can be cal- 
culated by integrating C/T, where C is the heat capac- 
ity. We neglect the distinction between constant volume 
and constant pressure for low compressibility solids and 
assume that the dominant entropy contribution comes 
from the lattice phonons. Note that at low-T the Debye 
approximation C = (127r 4 /5)ks(T/0) 3 becomes exact, 
where 9 is the Debye temperature. Neglecting any tem- 
perature variation of 9, the entropy also varies as (T/9) 3 . 
A value of $or=292K has been reported experimen- 
tally [20] ■ As we have not located an experimental value 

1/2 

for 9 c t we estimate 9 c t ~ (-BctV-Bor) #or=340K. 
Note that A(6»- 3 ) « (-3A6»)/6» 4 . 

Approximating AV and A9 as constants allows us to 
integrate the Clausius-Clapeyron equation, yielding 



P C (T) = P c ° - 



3tt 4 



k B A9\ (T 
AV 



(2) 



Putting in our numerical values, ksA9/ AV—Q.&Q GPa, 
and P c °=0.36 GPa, we can invert Eq.[2]to find the T C (P) 
phase boundary: 



T C (P) = 117 K x (P- 0.36) 



1/4 



(3) 



in units of Kelvin. Eq. [3] implies that the boundary of 
the collapsed tetragonal phase rises vertically from its low 
temperature limit of 0.36 GPa, then bends sharply to the 
right towards higher pressures, in qualitative agreement 
with experiment @, 0] • The boundary curves to the right 
towards higher pressure because the higher bulk modu- 
lus of the cT phase reduces its entropy and hence raises 
its Gibbs free energy, G=E+PV-TS, relative to the OR 
phase. At high temperature, it thus requires a higher 
pressure to favor the lower volume, but lower entropy, 
collapsed tetragonal phase. 

Experimentally the collapsed tetragonal cT phase bor- 
ders the orthorhombic OR phase at low temperatures, 
but it borders the non-collapsed tetragonal T phase at 
high temperatures. In fact, the high temperature T phase 
resembles the OR phase in terms of its lattice parameters 
and even exhibits magnetic moment fluctuations match- 



ing the striped antiferromagnetic structure with correla- 
tion lengths of order 6-8 A [8j . The T phase is essentially 
the OR phase with a loss of long-range order in magne- 
tization leading to a loss of orthorhombicity (i.e. a = b), 
and hence the extrapolation of the cT-OR boundary to 
higher temperature and pressure should remain a good 
first approximation to the cT-T phase boundary Q . 

We thank A. Goldman for useful discussions and ac- 
knowledge the support of Aspen Center for Physics, 
where part of the work was carried out. 
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